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ABSTRACT 

Context. Early-type dwarf galaxies are the most common type of galaxies observed in the Universe. Their study has important 
cosmological implications because according with hierarchical galaxy evolution theories they are the progenitors of brighter galaxies. 
Nevertheless, the origin of this kind of systems is still not well understood. 

Aims. The aim of the present work is to investigate whether the different locations of dwarf galaxies with respect to ellipticals in the 
face-on view of the fundamental plane could be due to the transformation of bright disc galaxies in low-mass systems by harassment. 
Methods. We have run high-resolution N-body numerical simulations to test the tidal stripping scenario of the dE galaxies. The 
present simulations modelled several individual tidal stripping events on initial disc-like galaxy models with different bulge-to-disc 
mass ratios. 

Results. The models have shown that tidal stripping is a very efficient mechanism for removing stars and dark matter particles from 
galaxies, specially from their outer parts. The particles of the disc and halo components were easily stripped, while the bulge not. 
Thus, the scale length of the discs were 40-50% shorter than the initial ones. Prograde tidal interactions create tidal features like stable 
bars in the discs of the galaxies. In contrast, bars are inhibited in retrograde encounters. After several tidal interactions the galaxy 
remnants looks like a dwarf spheroidal system. The final position of the low-mass systems in the face-on view of the fundamental 
plane depends on the initial conditions of the simulations. Thus, simulated galaxies with initial large B/D ratios are closer to the 
face-on view of the fundamental plane defined by bright E and bulges of early-type galaxies. Nevertheless, galaxies with initial small 
B/D ratio are located, after four fast tidal encounters, at the position of dE galaxies in the face-on view of the fundamental plane.The 
final position of the remnants in the FP do not depend on the orbit configuration of the encounters. 

Conclusions. We conclude that fast galaxy-galaxy interactions are efficient mechanisms for transforming bright galaxies in dwarf 
ones. Indeed, the different location observed between Es and dEs in the face-on view of the fundamental plane can be explained by 
the formation of dwarf galaxies by harassment of late-type bright ones. 

Key words. Galaxies: evolution - Galaxiesxlusters - Galaxies:interactions 
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1. Introduction 

Early-type dwarf galaxies are low luminous (Mb > -18) and red 



systems mostly located in galaxy clusters ( |Phillipps et al. 1998 
Prac y et al. 2004[ ). The study of these objects has cosmological 
interest because according with the hierarchical theory of galaxy 
evolution, dwarfs are the building blocks of the brightest galax- 
ies ([Kauffmann & Whit e 19931 ). Nevertheless, their origin is still 
matter of debate. Three are the most important formation sce- 
narios of early-type dwarf galaxies: gravitational collapse simi- 
lar to bright elliptical (E) galaxies, ram-pressure stripping from 
gas-rich low mass irregular galaxies, and tidal debris left from 
interactions between galaxies or with the gravitational potential 
in galaxy clusters. 

In the standard picture, early-type dwarf galaxies are 
formed from the gravitational collapse of primordial density 
fluctuations, similar to bright E galaxies. Once the first stars 
are formed, mechanisms of energy feedback into the interstellar 
medium have been proposed in order to regulate the subsequent 
star formation or even change the structure of the galaxies 
dDekel & Silk 1 986 ). This standard picture is mainly supported 
due to several tight correlations between the structural param- 
eters of bright and dwarf early-type galaxies. Those relations, 
discovered in Virgo, Fornax and Coma clusters, include those 



between colour and magnitude (Sandage & Visvanat han 1978|l, 
radius and luminosity jGraham & Guzm an 2003 ; 
Gutierrez et al. 200 4} |Aguerri et al. 2004}, surface brightness 
and luminosity (Kormendy 1977| |Binggeli & Cameron 1991 



Aguerri et al. 2005a IGraham & Guz man 2003J7 or Sersic 
shape parameter and luminosity ( Grah am & Guzm an 2003 : 



Aguerri et al. 2005a I. These tight correlations could indicate a 



common origin between bright E and early-type dwarf galaxies. 
However, those common relations have also been explained 
by some authors as being more due to the uniformity of their 
dark matter halos than to a common formation mechanism 
dNavarro et al. 19971 >. 

The second formation scenario is based on the presence 
of common stellar distributions of early-type dwarf and dwarf 



irregular (dlrr) galaxies dLin & F aber 198 31 |Korme ndy 1985 
ICaldwell & Bothun 19871 |Binggeli & Cameron 1991 

Ivan Zee et al. 2004al l. The role played by the environment in 
this formation scenario is crucial. Thus, ram pressure stripping 
can disrupt the gas from gas rich low-mass galaxies and trans- 
form them into early-type dwarf galaxies (Bo selli et al. 2 008: 
Sanchez-Janss en et al. 20 08 ). Repeated tidal shocks suffered by 
a dwarf satellite galaxy can also remove the kinematic signature 
of a dlrr galaxy into an early-type dwarf (Mayer et al. 2001a; 
Mayer et al. 2001b ). This formation mechanism has been rein- 
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forced recently by kinematics studies of early-type dwarf galax- 
ies in the Virgo cluster dPedraz et al. 20021 IGeha et al~2 002 ; 
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Geh a et al. 20031 Ivan Zee et al. 2004bl ICorsiniet al. 20071 ). 
It was found that most of the early-type dwarf galaxies are 
rotationally supported similar to dlrrs. But, a scenario where 
all early-type dwarf galaxies evolve from dlrr via disc fading 
do not however seem possible, because many of early-type 
dwarfs in the Virgo and Fornax clusters are brighter than the 
dlrr (Conselice et al. 20011 1. Another important point against 
this scenario is the increase of the metallicity measured in the 
old stellar populations in early-type dwarf galaxies with respect 
to dwarf irregulars of similar luminosity (Gre bel et al. 2003] ) . 
Although this objection was recently explained by Boselli et al. 
(2008). 

In the harassment scenario the cumulative effects of tidal 
fast encounters between galaxies and with the gravitational po- 
tential of the galaxy cluster can produce a dramatic morpho- 
logical transformation from bright spiral galaxies to dwarfs 
(Moo re et al. 19961 pastropietro et al. 2005| >. The effects of the 
harassment can be seen in the scale length of the discs of bright 
spiral galaxies. Thus, Aguerri et al. (2004) and Gutierrez et 
al. (2004) found that the scale length of the discs of bright 
spiral galaxies of the Coma cluster was shorter than simi- 
lar galaxies in the field. The disrupted stars could form part 
of the so-called diffuse light observed in some nearby clus- 
ters such as Virgo (Arnaboldi et al. 2002; Feldmeier et al. 2004; 



Aguerri et al. 2005b[ IWilliam s et al. 2007 ), Coma ( Gerhard et 
al. 2005, 2007) and Fornax dTheuns & Warren 19971 l . However, 
the major part of this diffuse cluster component seems to be as- 
sociated with the merger history of the brightest cluster galaxies 
rather than with disrupted dwarf galaxies (Mur ante et al. 20071) . 
Nevertheless, some bright spiral galaxies can lose part of their 
stars by fast interactions with the cluster potential and transform 
into a hot dwarf-like object similar to the observed early-type 
dwarf galaxies (Mastropietro et al. 2005 1. As previous dwarf 
formation scenario has also some objections. In particular, one 
prediction of the harassment model is that dwarf galaxies should 
lie in low surface brightness tidal streams or arcs. However, 
Davies et al. (2008) have found no evidence for any such fea- 
tures in a sample of dwarf galaxies in the Virgo cluster. 

In a recent series of papers, Lisker et al. (2006a; 2006b; 
2007) have proposed that the cluster dwarf galaxies population 
is composed of different sub-categories of objects, of which not 
all have been formed in the same way. The different formation 
mechanisms among the population of early-type dwarf galaxies 
could be seen in different positions on the fundamental plane 
(FP). The analysis of the location in the FP of harassed bright 
spiral galaxies is one piece of information lost in this puzzle, 
and will be analysed in the present paper. 

The FP is a tight relation between the surface brightness, 
effective radius and central velocity dispersion followed by 
hot systems such as elliptical galaxies and bulges of spi- 
rals. One of the projections of this plane is the so-called 
Faber-Jackson (IFaber & Jackson 19761 FJ) relation. It has 
been observed that dE and bright E galaxies have dif- 
ferent FP Bender et al. 19921 IPeterson & Caldwell 19931 
and Faber-Jackson relations (IHeld et al. 19921 |Tonry 1981 ; 
IDavies et al. 19831 IMatkovic & Guzman 20051 ) . In particular, 
the face-on view of the FP shows that dE are located in a 
disjoint area from bright E galaxies. It was also observed 
that bright E galaxies follow the FJ relation which can be 
expressed as L oc cr", where n » 4 (IFaber & Jackson 19761 
Sarg ent et al. \911\ ISchechter & Gunn 19791 [S chechter 19801 
Tonry 1 98 1 [ ITerlevich et al. 1 98 1 1 IForbes & Ponman 1 999 ), 
while dE galaxie s show n » 2 (|Tonry 198 1 1 IDavies et aT7 l983; 
IHeld et al. 19921 IMatkovic & Guzma n2005l These observa- 



tional evidences suggest a different origin and/or evolution of 
these two kinds of systems. Thus, Matkovic & Guzman (2005) 
demonstrated that the different FJ relation observed between 
bright E and dE can not be explained with the presence of 
rotation in the dE sample. They conclude that faint early-type 
galaxies in Coma had to be formed at least 6 Gyr ago and over a 
short 1 Gyr period. The different FP displayed by bright E and 
dE was also investigated by de Rijcke et al. (2005). They pointed 
out that internal processes could produce those differences. 
They concluded that semi-analytical models that incorporate 
quiescent star formation combined with post-merger starbursts 
and the dynamical response after supernova-driven gas-loss, 
were able to reproduce the position of dE in the FP. These 
two previous works indicate that the differences in the FP and 
FJ relation could be explained by different internal evolution 
between bright E and dE galaxies. Nevertheless, few attention 
has been given until now to the analysis of these differences due 
to other origins. 

The origin and evolution of the FP of elliptical galax- 
ies has received considerable attention from N-body modellers 
HCapelato et al. 1995| IBekki 19981 |Aceves & Velazquez 2005 1. 
Gonzalez-Garcfa & van Albada (2003) and Nipoti et al. (2003) 
showed that the FP is not severely affected by a few merging 
events. This results have been confirmed by several latter stud- 
ies (e.g. Boylan-Kolchin, Ma & Quataert 2006) including self- 
consistent cosmological simulations (Onorbe et al. 2006, 2007). 
However the location of the different galaxy types, like dEs, in 
the FP has received much less interest. 

The aim of the present work is to investigate whether the 
position of the dE in the FP can be explained by the harassment 
scenario, which posits that dEs stem from perturbed and stripped 
late-type disc galaxies that entered in clusters and groups of 
galaxies some Gyr ago. The structure of the paper is as fol- 
low: the sample of galaxies used and the FP of bright and low- 
luminous galaxies are shown in Section 2; N-body numerical 
simulations of harassed bright galaxies are presented in Section 
3; our discussion and conclusions are given in Sections 4 and 5, 
respectively. 

2. The Fundamental Plane 



It is well known that dynamically hot systems (elliptical galaxies 
and bulges of spiral galaxies) are not scattered randomly in the 
three-dimensional space spanned by the luminosity, half-light ra- 
dius, and velocity dispersion (Dressier et al. 1987; Djorgovski & 
Davis 1987; Bender, Burstein & Faber 1992; Falcon-Barroso, 
Peletier & Balcells 2002; Bernardi et al. 2003) . Those objects 
are located in a plane called the Fundamental Plane given by: 



log(r e ) = Qflog(cr ) +(3(fl e ) + J 



(1) 



where a, /3, and y are constants, and r e , (jJ. e ), and <x are the 
half-light effective radius, mean surface brightness and central 
velocity dispersion of the galaxies, respectively. 

The FP can be directly derived from the virial theorem: 



<v 2 > = GM/r, 



(2) 



where M is the total mass of the system, r g is known as the grav- 
itational radius and (v 2 ) is the mean square velocity of the stars 
in the system. To relate these quantities to observational ones, 
Capelato et al. (1995) introduced: 



C v = {v 2 )/cr 2 



C r — Fgl ?~e 



(3) 
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Fig. 1. Fundamental plane in the k representation of the differ- 
ent types of galaxies in the Coma cluster: bright E (empty cir- 
cles), dE (black circles), early-type bulges (red squares), late- 
type bulges (blue triangles), and dSO bulges (red circles). The 
full lines shows the fits of the FP in the /e-space defined only by 
bright Es. 
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Fig. 2. Residuals of the face-on fundamental plane in the k rep- 
resentation of the different types of galaxies in the Coma cluster: 
bright E (empty circles), dE (black circles), early-type bulges 
(red squares), late-type bulges (blue triangles), and dSO bulges 
(red circles). 

Introducing these parameters in eq (1), and using (jx e ) = 
-2.5 log(^r), being L the luminosity of the object, the FP can 
be written as: 



Thus, the virial theorem predicts that a — 2 and (3 = 0.4, 
while the observed values are a = 1.24 and /? = 0.328 
(J0rgensen, Franx & Kjaergaard, 1996). Equations (1) and (4) de- 
termine two oblique planes defining an angle, the so-called tilt 
of the FP. From eq. (4) can be seen that the tilt of the FP is due 
to differences in the structure or in the M/L ratio of the galaxies. 

Another formulation of the FP is the so-called K-space 
(Bender et al. 1992) expressed by: 

1 , 1 ho-l 1 crl 

ki = -= log(r e c^); k 2 = — log(^); * 3 = -= log(-2-) (5) 
V2 V6 r e V3 hr e 

Using these definitions and the virial theorem, k\ is expected 
to be a measure for the mass of the galaxy, k 2 is mostly sensitive 
to surface brightness, and depends on the mass-to-light ratio. 
In this parameterisation, the k\ - kt, plane is an edge-on view of 
the FP, while the k\ - k% plane gives the face-on view of the FP. 

Most of the studies about the FP have been centred on bright 
early-type galaxies in nearby galaxy clusters (see D'onofrio et al. 
2006, and references therein). Notwithstanding, the location of 
dwarf galaxies in the fundamental plane has not been extensively 
analysed after the early works by Nieto et al. (1990), Bender et 
al. (1992), and Guzman et al. (1993). However, Bender et al. 
(1992) and Guzman et al. (1993), based on a small sample of 
dwarf galaxies, found that dwarf galaxies were located in a dif- 
ferent position in the face-on view of the FP defined in the k- 
space. The distribution of the dwarf galaxies were perpendicular 
to the relationship defined by luminous elliptical galaxies in the 
K\ - K2 plane. This trend was also recently confirmed on a larger 
sample by De Rijcke et al. (2005). 

In the present work we have analysed the fundamental plane 
of a large sample of bright and dwarf galaxies located in the 
Coma cluster. The sample is formed by 215 galaxies of differ- 
ent types: 35 bright E, 113 dE, 12 dSO, 23 bulges of early-type 
spirals, and 1 1 bulges of late-type galaxiefl The photometric 
parameters of the galaxies (r e and [i e ) were obtained from the 
surface brightness decompositions on these objects carried out 
by Aguerri et al. (2004, 2005a), and the central velocity disper- 
sion of the galaxies (<x ) was taken from Matkovic & Guzman 
(2005) and Smith et al. (2004). 

Figure 1 shows the edge-on and face-on representations of 
the FP in the /c-space for the galaxies in the Coma cluster. 
We have fitted, in both orientations, the FP determined only 
by bright E galaxies. The scatter around the edge-on view of 
this FP depends on the galaxy type, being: 0.07, 0.11, 0.11, 
0.07, 0.23 for bright E, dE, dSO, early-type bulges and late-type 
bulges, respectively. Notice that low-luminous galaxies and late- 
type bulges have larger dispersion than bright E and early-type 
bulges. This is in agreement with previous findings by Nieto 
et al. (1990). The segregation between low-mass systems and 
bright E galaxies is clear in the k\ - k 2 plane (face-on view of the 
FP; see Fig. 1). As it was observed in previous works (Bender 
et al. 1992; Guzman et al. 1993; De Rijcke et al. 2005), dwarf 
galaxies and late-type bulges have different location on the face- 
on view of the FP than bright E and early-type bulges. This dif- 
ference is more clear in Fig. 2, which shows the residual of each 
galaxy respect to the fitted face-on view of the FP defined by 
bright E galaxies. In Fig. 2 we have scaled the symbols of the 
dwarf galaxies by their cr value. Notice that dwarf galaxies with 
low values of cr are the galaxies located further away from the 
FP defined by bright E. 



log(r e ) = 2 log(o- ) + 0.4<^ e > + log(C„C r ) - log(-) + y (4) 



1 See Aguerri et al. (2004, 2005a) for details on the galaxy classifi- 
cation 
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Fig. 3. Circular velocity curves for the initial conditions of the 
different simulations. 

Table 1. Input number of particles and mass fraction of each 
component of the different models. Particle number is given in 
K-particles (xlO 3 ). 



Model 


# 


$halo 


$disc 


$ bulge 


*>/„,„ 


Mji, c 


TID1 


370 


300 


60 


10 


1/8 


1/12 


TID2 


540 


275 


175 


90 


1/7.5 


1/5 


TID3 


540 


275 


175 


90 


1/5 


1/2 


TID4 


540 


275 


175 


90 


1/5 


1/2 


TID5 


400 


310 


90 




1/6 




TID6 


400 


310 


90 




1/6 




TID7 


400 


310 


90 




1/6 




TID8 


400 


310 


90 




1/6 




GST1 


370 


300 


60 


10 


1/8 


1/12 


GST2 


540 


275 


175 


90 


1/7.5 


1/5 



The aim of the present work is to understand if the segrega- 
tion between bright and dwarf galaxies in the face-on view of 
the FP could be explained by the transformation of bright disc 
galaxies in low-mass systems by fast galaxy-galaxy encounters. 



3. Numerical Simulations 

We have performed several N-body simulations to test the tidal 
stripping formation scenario for dEs. The present simulations 
are not cosmological ones. They model several individual tidal 
stripping events happening on different initial disc-like galaxy 
models. The simulations are pure N-body models with no gas 
content. Nevertheless, we have also simulated the effect of gas 
stripping on two models by removing 10% of the luminous par- 
ticles at the beginning of the simulations. After each tidal in- 
teraction or gas stripping event we allow to the system to reach 
virial equilibrium and then we have investigated its position in 
the K-space. 

3.1. Initial conditions and units of the simulations 

The initial galaxies of our simulations consisted on disc galaxies 
generated by the Kuijken & Dubinsky (1995) disc-bulge-halo 
algorithm. The galaxy has a King (1966) model for the bulge 
(if present), a Shu (1969) for the disc and an Evans (Kuijken & 



Dubinsky 1994) model for the halo, being all components mod- 
elled live. TableQ]gives the details of the different initial models 
of the simulations. 

We have simulated galaxies with different bulge-to-disc 
(B/D) mass ratios in order to test different starting points for 
the simulations. Thus, the initial B/D mass ratios were 1:12 
(TID1 and GST1 models), 1:5 (TID2 and GST2 models), and 1:2 
(TID3 and TID4 models). In addition, four more models (TID5, 
TID6, TID7 and TID8) were pure disc galaxies, where no bulge 
was considered in the initial conditions. Assuming the well- 
know observational correlation between morphological type and 
B/D (e.g. rSimien & de Vaucouleurs 19861 | Aguerri et al. 2004| 
IMendez-Abreu et al. 20081 we can say that TID3 and TID4 
could represent an early-type Sa galaxy; TID2 and GST2 show 
an Sb galaxy; TID1 and GST1 are similar to Sc galaxies; and 
TID5-TID8 can be compared with an Sd/Sm galaxy. 

We used model units to run the simulations. This units take 
the total mass of the initial galaxy as unity, the disc scale length 
as 0.2, and the Constant of Gravity G=l. However, in order to 
compare with real galaxies we need to convert the model units 
to physical ones. This was made by defining a set of conversion 
units to compare our simulations with observations. For bright 
galaxies, the scale length of the disc is independent of the mor- 
phological type, being i » 3 kpc (see e.g. Aguerri et al. 2004 
and references therein). This means that we can take a conver- 
sion factor of [L] = 15 kpc for all simulations as the unit of 
length. In contrast, different galaxy types show different rota- 
tion curves. Thus, the maximum velocity (V max ) of Sa, Sb, Sc, 
and Sd galaxies are ~300, ~250, ~180 and ~150 km/s, respec- 
tively (Roberts 1978; Rubin et al. 1985). This means that the 
velocity conversion factor for our simulations were [V] = 496 
km/s for TID3 and TID4; [V] = 345 km/s for TID2 and GST2; 
[V] = 230 km/s for TID1 and GST1; and [V] = 175 km/s 
for TID5-TID8. Figure 3 shows the initial rotation curves of 
the different simulations. Notice the different V max and rotation 
curve central gradients of the simulations. Thus, TID5-TID8 
models shows the smaller V max and the less steep central gra- 
dient in the rotation curve as it is observed in late-type galaxies. 
In contrast, TID3 and TID4 models show the largest V max and 
steeper rotation curve as correspond to Sa galaxies with large 
central mass concentration (large bulges). Assuming a time unit 
[T] = 2.4 x 10 7 yr, the total mass of the galaxies turned to be: 
[M] = 8.58 xlO n M for TID3 and TID4; [M] =4.15xl0 u M o 
forTID2 andGST2; [M] = 1.84 x 10 n M Q forTIDl andGSTl; 
and [M] = 1.07 x 10 n M Q for TID5-TID8. Thus, in agree- 
ment with observations, the more massive galaxies correspond 
to those models simulating early-type galaxies. 

3.2. Galaxy-Galaxy encounters 

The aim of the present simulations is to model fast encounters 
of our initial disc galaxy with another system of similar mass. 
Those galactic interactions would be similar to the fast galaxy- 
galaxy encounters predicted in high density environments as 
galaxy clusters (Moore et al. 1996). This kind of encounters 
has the peculiarity that the effective duration of the encounter 
is shorter than the characteristic crossing time of the major- 
ity of their stars. Moreover, the impact parameter of our fast 
encounters will be much larger than the typical galaxy scales. 
Therefore, this kind of galaxy encounters can be modelled using 
the impulse tidal approximation (see Binney & Tremaine 1987, 
§7. 2(d)). Thus, we can assume that the galaxies have no time to 
reorganise during the encounter, hence the energy injection is all 
kinetic; reorganisation occurs after the encounter has finished. 




Fig. 4. Evolution of the tidal events in TID1. Rows give the three views of the system after the different tidal harassment events. 
Colours indicate the different constituents of the initial galaxy: black-dark halo, green-disc and red-bulge 
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Fig. 5. Evolution of the tidal events in TID3. Rows give the three views of the system after the different tidal harassment events. 
Colours indicate the different constituents of the initial galaxy: black-dark halo, green-disc and red-bulge 
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The commonly used formula for the tidal impulse approx- 
imation does not take into account the extended nature of both 
bodies. Nevertheless, our encounters could not be considered be- 
tween point-like galaxies. To take this into account, we have ap- 
plied the formalism proposed by Gnedin, Hernquist & Ostriker 
(1999) giving that the variation in the velocity (AV2), angular 
momentum (A/2) and energy (AE2) of the particles in the mod- 
els are: 



Table 2. Half mass radius for each of the components for the 
initial and final models of our simulations. 



GM] Rperi 1 



AE 2 (R) 



AJi(R) 



b 2 V R„ 



G 2 M? , 
' b 4 V 2 

GMi 



(6) 

(7) 
(8) 



where Mi is the total mass of the perturber, b is the impact pa- 
rameter of the orbit, R per i is the pericenter distance, V is the rela- 
tive velocity at pericenter, R ma x is the cut-off imposed in our ini- 
tial model for the perturber galaxy, and I22 is the smallest eigen- 
value of the inertia tensor. A detailed discussion of these for- 
mulae and their implications can be found in Gonzalez-Garcfa, 
Aguerri & Balcells (2005). 

The fast encounters were simulated injecting energy, accord- 
ing with the previous expressions, to each particle of the mod- 
elled galaxy. Stars and dark matter particles in the modelled 
galaxy absorb energy (AE =s 0.5(AV2) 2 ) as well as angular mo- 
mentum and may become unbound after the interaction. Since 
the velocity and energy change depend on the distance to the 
central parts, the particles with higher binding energy (E^md) are 
the ones that are going to be easily stripped away. After each 
tidal interaction we allow the system to reach virialization in the 
innermost regions by evolving the galaxy in isolation for more 
than 25 time units. This amounts for more than 50 half-mass ra- 
dius crossing times of the bulge. For the different runs we have 
taken b = V2, V = 2, R peri = R e and M\ equal to the total 
mass of the initial galaxy at each tidal interaction, so we model 
equal-mass encounters. A total number of 4 tidal interactions 
were simulated for each model. 

We have run 10 different simulations (see Tab. 1). The mod- 
els called TID1, TID2, TID3 and TID5 differ on the bulge-to- 
disc ratio (B/D). Thus, they have B/D = 1/12, 1/5, 1/2, and 
respectively. The relative orientation of the initial disc with 
respect to the orbital angular momentum was (6, </>) = (0,0). 
Models TID4 and TID6 were replicas of TID3 and TID5, chang- 
ing the orientation of the initial discs with respect to the assumed 
orbit, now being (60, 30). Models TID7 and TID8 are also iden- 
tical to model TID5 but with retrograde orbits with respect to the 
disc spin; TID7 has (6, (p) = (180, 0), while TID8 has a less sym- 
metric configuration (9, </>) = (135, 0). The tidal features of these 
systems were somewhat different from TID3 and TID5, but the 
amount of material stripped and the general structure of the rem- 
nants were basically the same. We have also run two more mod- 
els (GST1 and GST2) simulating an initial gas stripping. Those 
simulations have identical initial conditions to TID1 and TID2, 
respectively. The initial gas stripping was simulated by randomly 
removing 10% of the disc particles of the models. Gas-stripping 
depends also on the distance from the centre of the galaxy (see 
the hydrodynamical simulations by Quilis et al. 2000), but here 
we adopt an extreme situation where at a single go all gas in the 
galaxy is stripped away. After they reach virialization we applied 
the first tidal encounter. 
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Fig. 6. Change of mass in each tidal harassment event. Upper 
panektotal mass, lower panel: luminous mass. 



Simulations were run using the GADGET2 (Springel 2005) 
hydrodynamical code, in the tree-code configuration, on the 
Beowulf cluster of the Instituto de Astroffsica de Canarias (IAC), 
the Mare Nostrum supercomputer and the La Palma node of 
the National Supercomputing Center (BSC). A typical run with 
5 x 10 5 particles took of the order of 5 x 10 4 s in 16 CPUs. 
We used different softening parameters for the different parti- 
cles, typically e ~ 1/10 of the half-mass radius of each system. 
An opening parameter of 6 = 0.8 and quadrupole terms were 
included in the force calculation. 



3.3. Results 

Figures |4]and |5]present the initial and evolved states after each 
tidal interaction for the TID1 and TID3 simulations, respec- 
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tively. These two models represent two extreme cases because 
they have the largest (TID3) and smallest (TID1) initial bulges. 
In both cases, tidal stripping removes first the outer parts of the 
halo, after two tidal interactions a fair amount of disc particles 
are also expelled. Models with large bulges (TID2, TID3) show 
that almost the whole disc was removed after 4 tidal interac- 
tions. These models have at the end a disc-to-bulge ratio close to 
1:1, leaving the small number of still bounded disc particles in a 
roundish component (see Fig . O . In contrast, the model with the 
smallest initial bulge (TID1) still presents a short and thick disc 
after 4 tidal encounters (see Fig.[4]i. These truncations have been 
quantified measuring the half mass radius of the bulge (r e j,ui ge ), 
disc {r e c ii SC ) and halo (r e haio) as showed in Table 2. Notice that 
while the disc and halo are strongly truncated in all simulations, 
the bulge enlarge its half mass radius specially for the models 
with the largest B/D ratio (TID3 and TID4). 



Table 3. Maximun rotation velocity and velocity dispersion of 
the initial and final models and mean VI cr inside the effective 
radius for the final model. 



Vmaxjnitial Vmax,fina\ @" maxjnitia! ^maxjinal C^/^")/ina/ 

(km s~') (kms~') (kms -1 ) (kms -1 ) 
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3.3.1 . Disc and Halo truncations 



The truncation of the disc and halo component results in an im- 
portant mass loss of the galaxies. Figure[6]presents the evolution 
of the total content of bound mass in these systems. We find 
that harassment is highly efficient in removing a large amount of 
the initial mass in all our simulations. The different encounters 
removes between ~ 50% to 80% of the total initial mass, and 
between 30% and 60% of the luminous mass. 

Figure [7] shows the binding energy vs. the radial distance to 
the galaxy centre after each tidal encounter for TID1 and TID3 
models. It is clear that the outermost particles are stripped first 
after each tidal interaction because they have larger binding ener- 
gies. Thus, disc and halo particles were easily stripped away con- 
trary to bulge particles. In both models shown in Fig.|7]the sim- 
ulated galaxy-galaxy encounters do not provide enough energy 
to extract the bulge particles. But, some bulge particles occupy 
larger binding energies previously occupied by striped disc and 
halo particles producing an inflation of the bulge. This behaviour 
is common to all the simulations, being more prominent in those 
with large B/D-ratio. According with this, we can say that the 
bulge survive in those simulations with initial bulge. Moreover, 
at the end of the simulations, the hot component of galaxies with 
initial bulges is basically formed by the initial bulge of the galax- 
ies. 



3.3.2. Tidally induced bars 

The stripping remnants show a number of tidal features like ex- 
tended tidal tails, rings, etc, that indicate that the system has suf- 
fered such interactions. 

Some systems appear to present a bar-like structure in the 
disc (see Figure [8]). This bar-like structures induced by tidal in- 
teractions are formed in the first tidal encounter and are stable 
until the end of the simulation. The bars are more prominent in 
those models with lower central mass concentration, i.e. smaller 
bulge component, and a prograde rotation. The bar formation de- 
pends on the orbit configuration of the encounters. Thus, retro- 
grade orbits seem to inhibit bar formation. Only the models with 
the largest bulge component (TID3 and TID4) do not show any 
bar feature in the disc. This is in agreement with numerical simu- 
lations about bar formation and evolution which find that central 
mass concentrations inhibits bar formation in discs. (Pfenniger 
& Norman 1990; Shen & Sellwood 2004; Athanassoula et al. 
2005). 



3.3.3. Kinematical features of the remnants 

Figure [9] shows the rotation curve and the velocity disper- 
sion profile of the initial and final remnants of all simulations. 
Table 3 gives the maximum rotation and velocity dispersion 
of the initial and final models. Notice that the rotation veloc- 
ity of all models are significantly smaller in the external re- 
gions at the end of the simulations than in the initial condi- 
tions. In contrast, the external region of the galaxies show also 
in all simulations larger velocity dispersion than at the begin- 
ning. This is caused because tidal interactions are an efficient 
mechanism for heating the thin disc of spiral galaxies (see e.g. 
|Aguerri et al. 200 1| lEliche-Moral et al. 2006b . The central ve- 
locity dispersion at the end of the simulations depends on the 
models. Those models with an initial small bulge show very 
similar central velocity dispersion at the beginning and at the 
end of the simulations (TID1, TID2). Models with large bulge 
at the initial conditions (TID3 and TID4) show smaller central 
velocity dispersion at the end. This is due to the final bulges of 
this models suffered an expansion. On the other hand, the mod- 
els without initial bulge (TID5-TID8) show larger central veloc- 
ity dispersion at the end, again due to the tidal heating of the 
disc and the absence of a stabilising central concentration like a 
bulge. Retrograde models seem to present slightly higher values 
of the central velocity dispersion. Table 3 also shows the V/cr 
ratio measured within an aperture equal to r e . Notice that the 
final remnants are in all cases dominated by random motions. 
Similar kinematics features have also been seen in other simula- 
tions of fast galaxy encounters in the literature (see Kronberger 
et al. 2006). 

The kinematical characteristics of the models with simulated 
tidal stripping (GST1 and GST2) do not differ largely from those 
similar simulations without tidal stripping (TID1 and TID2). 
This is because the gas-particles only account for a few per- 
cent (10% in our simulations) of the total mass of the disc in 
the galaxies. 



3.4. The Fundamental Plane of the tidal remnants 

We have analysed the location of the simulated galaxies in the 
FP after each tidal encounter. We have measured the required 
quantities ji e , r e and <x , in the simulations. 
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Fig. 7. Binding energy vs distance from the centre for the final snapshots of the simulations TID1 (up) and TID3 (bottom). The first 
panel in each case presents the state of the initial galaxy model. 
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Fig. 8. Simulation TID2, snapshot after the evolution of the second tidal harassment event in simulation TID2. Left and central panel 
show the luminous particles (bulge+disc). We observe the formation of tidal tails (left) and a prominent bar (centre). This bar can 
clearly be seen in the contour plot (right). 



3.4.1 . Decomposition of the surface mass density profiles 

Similar to real galaxies, we have performed the decomposition 
of the surface mass density profiles in the different galactic com- 
ponents of the initial modelled galaxy and the remnants after 
each interaction. This decomposition was done using one or two 
components depending if the systems are dominated by rota- 
tion or pressure in their most external regions. Thus, systems 
with v/cr < 1 at all radial distances were considered as hot sys- 
tems. In contrast, those systems with v/cr > 1 at large radial 
distances were considered as rotation supported systems in their 
outer most regions. Thus, the surface mass density profiles of 
the galaxies with a rotating component in their outermost re- 
gion were fitted with two components: bulge and disc. The bulge 



component dominates the surface mass density profiles in the in- 
nermost regions of the models, and was represented by a Sersic 
law dSersic 19681 I. This component represents the hottest part 
of the galaxy being dominated by random motions, and it will 
be called "bulge" similar to real galaxies. The disc component 
dominates the outermost regions of the mass density profiles of 
the galaxies, and was modelled by a Freeman exponential law 
(Freeman et al. 1970). The surface mass density profiles of those 
models dominated by random motions at all radial distances 
were fitted by one-component model described by a Sersic law. 

The free parameters of the models were fitted using a 
Levenberg-Marquard nonlinear fitting routine and were deter- 
mined by minimising the x 2 of the fit similar to real galaxies 
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Fig. 9. Rotation curves (left panels) and dispersion velocity profiles (right panels) for the initial (full lines) and remnants (dashed 
lines) of the simulations. 



(Trujillo et al. 2001, 2002; Aguerri & Trujillo 2002; Aguerri et 
al. 2004, 2005a). 

The decompositions of the surface mass density profiles pro- 
vided two of the parameters of the hot component of the galaxy 
needed for the FP, the effective radius (r e ) and the effective sur- 
face mass density (E e ). The third parameter (central velocity dis- 
persion) was measured from the velocity dispersion profiles of 
the galaxies in an aperture of radius equal to 0.2xr e of the hot 
galaxy component. All these FP parameters were obtained for 9 
different points of view in each system. 

As pointed out before, the discs of the galaxies are strongly 
truncated. This is reflected in the scale-length of the discs which 
is reduced in all cases up to 40-50%. In contrast, the variation 
in the effective radius of the bulges depends on the prominence 
of this structure in the initial model. Thus, the scale of bulges in 
models with small B/D ratio (TID1 and GST1) was very small. 
In contrast, galaxies with large B/D ratio (TID3 and TID4) have 
larger scales of the bulge at the end of the interactions. This be- 
haviours of the bulge and disc structural parameters produce that 
the B/D ratio of the tidally truncated galaxies is larger than the 
initial one. 

3.4.2. Location of tidal remnants in the FP 

To locate the modelled galaxies in the observed FP we need to 
scale them using an appropriate mass-to-light ratio. This M/L 



was obtained by imposing that the hot systems of the simulated 
initial galaxy and remnants lied on the edge-on view of the FP 
described by bright E galaxies. With this M/L we have compared 
the location of those systems in the face-on view of the FP. We 
have varied the imposed M/L according with the rms showed 
in the edge-on view of the FP of early-type bulges (TID2, TID3, 
TID4 and GST2) and late-type bulges (TID1 and GST1). The re- 
sult of this variation in the face on view of the FP is given as an 
error bar in the following FP figures. Figure [10] shows the loca- 
tion in the /e-space of the observed galaxies in the Coma cluster 
presented in § [2j together with the hot components correspond- 
ing to the TID1, TID2 and TID3 models. We find that those 
simulations with initial conditions similar to an early-type spi- 
ral (TID3 model) are located close to the face-on view of the FP 
described by bulges of early-type and E galaxies. In contrast, the 
remnants of the models with small B/D ratio at the beginning 
(TID1 model) are located in a different position in the face-on 
view of the FP defined by E galaxies, similar to dE ones. The 
remnants from TID2 are between TID1 and TID3 models. 

Models GST1 and GST2 have a gas-stripping event at the 
beginning of the simulations. This event is supposed to erase all 
gas from the galaxy in a rather short amount of time. This mod- 
ifies the profile of our initial conditions producing initial bulges 
with larger scales than in similar models without gas-stripping 
(TID1 and TID2 models). The final remnants for these models 
are located at different positions in the face-on view of the FP 
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than those from TID1 and TID2 simulations (see Fig. ITTb . In 
fact, we find a marginally better agreement with the locus of dEs 
in the model GST1. 

All these models with galaxies with bulge at the beginning 
can explain the location in the face-on view of the FP of dwarf 
galaxies with central velocity dispersion larger than w 50 km s _1 . 
But, dwarf galaxies with smaller central velocity dispersions are 
yet not well reproduced. 

The Ml L ratio for the four models without initial bulge com- 
ponent (TID5-TID8) were obtained in a different way. In this 
case the initial models represent pure disc galaxies, and it has 
no sense to locate these discs in the FP. We have obtained the 
M/L ratio by locating the initial discs of these models in the 
Tully-Fisher (Tully & Fisher (1977); TF) relation followed by 
spiral galaxies. We have the maximum rotational velocity for 
these discs, and following the Tully-Fisher calibration of Sakai 
et al. (2000) in the I band we have obtained the M/L. 

Figure[T2]shows the location of the remnants from TID5 and 
TID6 models in the edge-on (upper panel) and face-on (lower 
panel) for the FP. The final remnants from these models show 
smaller central velocity dispersions than the other simulations, 
and are located in the region occupied by dE with the smallest 
central velocity dispersions. These results do not depend on the 
orbit of the interactions. Thus, the models TID7 and TID8, fol- 
lowing retrograde orbits, show similar results. 



4. Discussion 

The simulations presented in the present paper indicate that 
fast tidal interactions are very efficient transforming spiral disc 
galaxies into dwarf like systems. This mechanism produces a 
large mass loss, specially from the outer parts of the galaxies, 
and locate the tidal remnants out from the face-on view of the 
FP described by bright E galaxies. Figure[T3]shows the residuals 
of the observed galaxies and the remnants of the present models 
respect to the face-on view of the FP described by E bright galax- 
ies. Notice the tidal remnants from the present simulations are 
located in the same position of the FP than dE galaxies. Galaxies 
with large B/D-ratios at the beginning are located after four tidal 
fast interaction near to the loci occupied by bulges of early-type 
galaxies. In contrast, galaxies with small B/D-ratio or without 
initial bulge are located at the end of the simulation at the po- 
sition of the dE galaxies with small central velocity dispersion. 
This indicates that dE galaxies can be formed by the tidal disrup- 
tion of late-type disc galaxies. In contrast, the tidal disruption of 
early-type spirals would produce a final galaxy similar to ellipti- 
cal or SO galaxies. 

This is not the unique mechanism that can explain the lo- 
cation of dwarf galaxies in the face-on view of the fundamen- 
tal plane. Semi-analytical models of dwarf galaxy evolution by 
supernova-driven-gas-loss can also explain the position of dwarf 
galaxies in the face-on view of the FP (see de Rijcke et al. 
2005). This could indicate that the origin of the dwarf popu- 
lation could be diverse. Thus, several past studies have shown 
that dwarf galaxies are not an homogeneous population of ob- 
jects. Recently, Lisker et al. (2007) have analysed a large sam- 
ple of dwarf elliptical galaxies in the Virgo cluster, conclud- 
ing that they can be grouped in three main families with dif- 
ferent dynamical and morphological properties. They found that 
non-nucleated dE show shapes like thick discs and are not cen- 
trally clustered. This indicates that they constitute a non-relaxed 
galaxy population, and could be formed by tidal interactions as 
our simulations show. But, how can we know the progenitors 



of dwarf galaxies?, are there any observational measurement for 
knowing the tidal origin of dwarf galaxies?. 



4. 1 . Bar parameters in dwarf galaxies 

As we mentioned in the previous sections, one of the most 
prominent features that tidal interactions induce in the discs of 
the galaxies is the presence of stable and prominent bars. Those 
bar features appeard only in those models with small bulge com- 
ponent and in interactions following prograde orbits. Retrograde 
orbits inhibits bar formation even in galaxies with small bulge. 
We can analyse the properties of these bars and compare them 
with similar properties of bars observed in bright galaxies with- 
out signs of interaction. Numerical simulations of bar formation 
find that the mean ratio between the bar length and the scale 
length of the galaxy discs (rj/A) is well constrained, being 1.0- 
1.2 (O'Neill & Dubinski 2003; Valenzuela & Klypin 2003). This 
is in agreement with observational measurements of this quantity 
in bright unperturbed barred galaxies (< rb/h >= 1.21 ± 0.08; 
Aguerri et al. 2005c). We have measured this ratio in the bars 
induced in our models by tidal interactions. Figure [14] shows 
the rb/h distribution on a sample of 14 bright barred galaxies 
(Aguerri et al. 2005c) and the values of this ratio for the bars 
formed in our simulations. It is clear than the tidally induced 
bars have larger rb/h values than observed bars. This could be an 
observational parameter that could be measured in dwarf barred 
galaxies in order to distinguish between tidally induced bars and 
those formed by internal discs instabilities. 

There are few measurements of the ratio rb/h in the literature 
for dwarf galaxies. Nevertheless, NGC4431 is one exception. 
This is a barred dwarf galaxy located in the Virgo cluster and 
analysed by Corsini et al. (2007). The surface brightness pro- 
file of this galaxy was decomposed in three components (bulge, 
disc and bar) resulting rb/h = 1.4, similar to the values obtained 
in our simulations. Therefore, the bar of this galaxy could be 
formed in a past tidal interaction. 



4.2. Structural parameters of tidal dwarf galaxies 

The analysis of the structural parameters of the dwarf galaxies 
can also give clues about their tidal origin. Photometric decom- 
position of bright and unperturbed galaxies have shown several 
correlations between bulge and disc structural parameters (e.g. 
de Jong 1996} IKhosroshahi et al. 20001 IMacArfhur et al. 20031 
Mollenhoff 2004, Hunt et al. 2 004; Me ndez-Abreu et al. 20081) . 



In particular, the ratio between the disc and bulge scales (r e /h) 
is correlated with the Sersic n parameter in real disc galaxies 
(e.g. MacArthur et al. 2003; Mendez-Abreu et al. 2008). Thus, 
those bulges with large n values also show large r e /h ratios. We 
have compared the correlation observed in real galaxies with our 
simulations. FigurefTBIshows the correlation between r e /h and n 
for a large sample of bright disc galaxies from the surveys of 
MacArthur et al. (2003) and Mendez-Abreu et al. (2008). We 
have also overploted in this figure the same quantities for our 
simulated dwarf-like objects. Notice that dwarf galaxies formed 
by fast tidal interactions are most of them outside the trend de- 
scribed by real galaxies. Thus, for a fixed value of n tidally 
formed dwarf galaxies have, in general, larger r e /h ratios than 
bright disc spirals. This is because, as we discuss in the previous 
section, the scales of the discs were strongly truncated due to 
tidal encounters. The search of outliers in the r e /h vs n relation 
could provide dwarf systems formed by transformed bright disc 
galaxies by fast tidal encounters. 
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We have the bulge and disc structural parameters of the 
dwarf barred galaxy NGC443 1 . This galaxy shows a Sersic in- 
dex n — 1.5 and r e /h = 0.44 (Corsini et al. 2007). According 
with Fig. [T3]this galaxy would be in the limit of the region de- 
scribed by the observed bright disc galaxies. This would be an- 
other indication that NGC443 1 could be a dwarf galaxy formed 
by stripped bright disc galaxies. 

4.3. Metallicity relations 

Early-type galaxies follow the well known luminosity- 
metallicity and velocity dispersion-metallicity relations (see e.g. 
Tremonti et al. 2004 and references therein). Thus, more metal- 
lic galaxies are those with larger central velocity dispersions 
and more luminous ones. Tidal stripping interactions should not 
change the metallicity of the galaxies. Indeed, we have also seen 
that the central velocity dispersion of tidal remnants are not very 
different from the initial conditions. In contrast, due to the large 
mass lost, tidally formed dwarf galaxies have lower luminosities 
than their progenitors. This would imply that low-mass systems 
formed by harassment of bright disc galaxies could be out of the 
luminosity-metallicity relation displayed by bright galaxies. In 
contrast, those dwarf galaxies could lie inside the relation be- 
tween the metallicity and their central velocity dispersion. 

4.4. Group and Cluster fast tidal interactions 

We have made some attempts to apply the tidal impulse approx- 
imation to model encounters with group and cluster potentials. 
Mastropietro et al. (2005) simulate the harassment of disc galax- 
ies falling onto a cluster environment by re-simulating a cluster 
formed in a cosmological simulation. We have tried to model 
the harassment of one of our bulgeless galaxies through the tidal 
impulse approximation. However, two important problems arise: 
The most typical orbital configuration from the cosmological 
simulations (Tormen 1997) are such that the injection of en- 
ergy is too low to produce significant effects in the discs. This 
is so because the typical orbits imply wide impact parameters on 
non-parabolic orbits. The second problem is that under this cir- 
cumstances no longer applies the impulse approximation (i.e., 
the interaction is not fast compared with the dynamic time of 
the galaxy), and the approach followed in this paper is no longer 
applicable. 

Nevertheless we have modelled one encounter with higher 
mass ratio (1:10). This model (called TID10) simulates the in- 
teraction of a bulgeless galaxy with the potential of a group. 
The orbital parameters chosen for this encounter were: b — 6, 
V = 4.6 and R per j = 1. This orbit is not a typical orbit (Tormen 
1997) but the parameters are such that the increment of energy 
to each particle is comparable to the ones in galaxy-galaxy in- 
teractions. In fact, the results are quite similar to the simulations 
TID5-TID8 (see figure 6 for the mass change). The final half 
mass radii of the halo and disc components were r ey haio = 0.532 
and r e ^isc = 0.370, respectively. The final maximum rotational 
velocity and central velocity dispersion were 37 km s _1 , and 
55.6 km s , respectively. All these values are similar to the ones 
shown in Tab. 2 for the models without bulge in the initial condi- 
tions (TID5-TID8). Indeed, the final position of these remnants 
in the FP was similar to similar models (TID5-TID8) showed 
in the present paper. However, we must stress that the orbital 
configuration is not typical at all, compared with those obtained 
from cosmological simulations (Tormen 1997). 



We would like also to point that the models showed in this 
subsection are very simple. The interaction of a galaxy with the 
global cluster or group potentials should take into account NFW- 
type tidal fields from the potential. This kind of interactions 
are out of the main scope of the present paper where galaxy- 
galxy interactions were shown. We are planning new cosmolog- 
ical simulations considering both the galaxy-galaxy and galaxy- 
cluster potential interactions. 

5. Conclusions 

In the present paper we have shown the FP of 215 hot systems 
(ellipcicals and bulges) located in the Coma clusters. The edge- 
on view of the FP shows that bright E galaxies and early-type 
bulges are located in the same plane showing similar scatter. In 
contrast, dwarf galaxies and late-type bulges have larger scat- 
ter. The differences are more pronounced in the face-on view of 
the FP. Thus, dwarf galaxies and bulges of late-type galaxies are 
located perpendicular to the relationship defined by E and early- 
type bulges. 

We have run high-resolution numerical simulations in order 
to explain the different location of dwarf galaxies in the face-on 
view of the FP. We have tested if fast galaxy-galaxy tidal inter- 
actions can transform bright spiral galaxies in early-type dwarf 
systems and can explain their different location in the FP. The 
simulations have disc-like galaxies with different bulge-to-disc 
mass rations as initial conditions. Thus, these initial disc galax- 
ies represent the full family of Hubble type galaxies, from early- 
to late-type ones. 

Tidal interactions turned out to be very efficient for remov- 
ing stars from the galaxies, specially from the outermost regions 
(disc and dark matter halo) where the stars are less bounded. 
In contrast, particles from bulges were not removed. After sev- 
eral tidal interactions the galaxies have lost 50-80% of their total 
mass and 30-60% of the luminous matter. The discs and haloes 
turned to be strongly truncated. In particular, the structure of the 
dies were largely affected. Their scale parameters were 40-50% 
shorter than the initial ones. Those models with the smallest B/D 
ratio retain a short and thick disc. In contrast, the model with 
the largest B/D ratio have no disc at the end, being a roundish 
dwarf system. The tidal interactions produce bar-like structures 
in the galaxy which are stable until the end of the models. The 
bars were more prominent in those models with smaller bulges 
at the beginning of the simulations. The formation of bar-like 
structures depends on the orbital configuration of the encounters. 
Thus, bars are easily formed in prograde interactions. In con- 
trast, retrograde orbits seem to inhibit the formation of bars even 
in those simulations without initial bulge component. Similar re- 
sults have been obtained for two models in which we have sim- 
ulated a gas stripping interaction before tidal encounters. 

We have analysed the location of the dwarf-like galaxies 
remnants in the face-on view of the FP. Their final positions 
depend on the initial conditions of the simulations. Thus, those 
dwarf galaxies with progenitors with large B/D ratio are closer 
to the location of early-type bulges and Es in the face-on view of 
the FP. In contrast, the remnants of late-type galaxies are located 
in the region occupied by dE in the face-on view of the FP. The 
final location of the remmants in the FP does not depend on the 
orbit orientation of the encounters. 

We have analysed several observational parameters in or- 
der to identify these tidal dwarf galaxies. The measurement of 
the rb/h ratio, the relation between r e /h and n, and the mass- 
metallicity and velocity dispersion-metallicity relations could 
provide information about the possible tidal origin of observed 
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dwarf galaxies. Therefore, tidal dwarf systems would have larger 
rb/h ratio than the observed one in bright barred galaxies. They 
would also be outliers in the r e jh-n and luminosity-metallicity 
relations. In contrast, they would follow the velocity dispersion- 
metallicity relation. 
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Fig. 11. Results for the models GST1 and GST2 presented in 
/f-space. Symbols and colours are as in Figur^3] 



Fig. 10. Results for the models TID1 (top), TID2 (middle) and 
TID3 (bottom) presented in the /(-space. The different types of 
galaxies in the Coma cluster are indicated as black small sym- 
bols: bright E (diamonds), dE (black circles), early-type bulges 
(red squares), late-type bulges (blue squares), and dSO bulges 
(red circles)). Colour symbols indicate the results for the mod- 
els. 
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Fig. 12. Results for the models TID5 and TID6 presented in the 
/f-space. Symbols and colours are as in Figur^3] 
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Fig. 14. Ratio between the length of the bars (rb) and the scale 
length of the discs (h) in observed bright unperturbed barred 
galaxies (histogram) and simulations with bars formed by tidal 
interactions (vertical lines). The different lines represent bars 
formed after 1 (full line), 2 (dotted line), 3 (dashed line) or 4 
(dashed-dot-dot line) tidal encounters. 



10.0 



c 1 .0 



••• •■ :•■ »• 

•• «• o- * 



0.0 



0.2 



0.4 

r./h 



0.6 



0.8 



Fig. 15. r e jh vs h for observed bright undisturbed galaxies (black 
points) and simulated dwarf galaxies formed by harassment of 
bright disc ones (diamonds). 



Fig. 13. Residuals of the face-on FP in the k representation of the 
different types of galaxies in the Coma cluster: bright E(empty 
circles), dE (black circles), early-type bulges (red squares), late- 
type bulges (blue triangles), and dSO bulges (red circles). The 
empty squares represent the simulated galaxies from models 
TID1-TID8, GST1 and GST2. 



